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An algebraic method is introduced for an analytical solution of the eigenvalue problem of the 
Tavis-Cummings (TC) Hamiltonian, based on polynomially deformed su(2), i.e. sii„(2), algebras. 
In this method the eigenvalue problem is solved in terms of a specific perturbation theory, developed 
here up to third order. Generalization to the A/'-atom case of the Rabi frequency and dressed states 
is also provided. A remarkable enhancement of spontaneous emission of M atoms in a resonator is 
found to result from collective effects. 
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I. INTRODUCTION 

We consider here the collective behavior of the system 
of M two-level atoms coupled to a single mode of the 
electromagnetic field in a resonator. The useful form for 
the atom-field interaction was proposed in the rotating 
wave approximation (RWA) by Tavis and Cummings . 
In their model M identical two-level atoms interact via 
dipole coupling with a single-mode quantized radiation 
field at resonance, so that the Hamiltonian is given by 



V = g{a'< ■ S^+a- S+) . 



(1) 



Here w is a frequency of the electromagnetic field and 
is the level splitting of the two-level atoms. The operators 
53, S± are collective spin variables of A/" two-level atoms. 
These operators are defined as 



E 



(2) 



where u's are Pauli matrices. They satisfy the su(2) alge- 
bra, a, are the annihilation and creation operators of 
the field. Due to historical reasons the Tavis-Cummings 
(TC) model is often called the Dicke model [||. We con- 
centrate here on the case of exact resonance, i.e. uj = loq. 
In this case the system exhibits a most interesting collec- 
tive behavior. For simplicity, time will be measured in 
units of the coupling constant g, i.e. we assume in the 
following that g — \. 



'Electronic address: 
t Electronic address: 
Electronic address: 
§ Electronic address: 



radeikoOmkk.ifmo.rii 




nirosli@nikk.ifmo.ri; 




:vl)in (.n:)li\"S.h'ii.fi 






ussi.timonenCSphys.jyu.fi 



The one-atom version of the TC model, i.e. the Jaynes- 
Cummings model(JC) has recently attracted much 
attention because of e.g. the spectacular experiments 
reported in jj]. In this work the properties of a dilute 
fiux of atoms excited to Rydberg states and interact- 
ing with a single quantum mode of the electromagnetic 
field were analyzed in a high-Q cavity. The results re- 
ported in ^] demonstrate that the physical assumptions 
underlying the JC model can indeed be experimentally 
realized. This triggered a lot of theoretical interest to- 
wards the one-atom micromaser (microlaser) (see ^ |^ 
and references therein). The analysis of the JC model is 
greatly simplified by the fact that the eigenvalue prob- 
lem of the J ay nes- Cummings Hamiltonian can be solved 
exactly. An algebraic approach to the solution of gen- 
eralized Jaynes-Cummings types of models can be found 
in ^, §). The approach to the Tavis-Cummings type 
models based on the Algebraic Bethe Ansatz is formu- 
lated in 0. 

An analytical solution of the A/'-atom case, i.e. the TC 
model, is much more complicated than that of JC model 
and is thus still far from completion. In this article we 
analyze this model and develop an analytical technique 
that is applicable to a variety of problems of TC type. 

To begin with we notice that Hamiltonian (^ belongs 
to a class of operators that can be expressed in the form: 



H = f{Ao) + g-{A+ + A^ 



(3) 



Here /(x) is an analytic function of x with real coeffi- 
cients, while the operators A± satisfy the commutation 
relations 



[Aa,A±] = ±A 



(4) 



We also assume the operator H to be self-adjoint. This 
means that, for any irreducible representation (irrep), the 
operators Aa,A± must satisfy the conditions 



(Ao)^ = Ao, {A J 



A 



+ ■ 



(5) 
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It is worth mentioning that Hamihonians of type (||) are 
also used in description of Raman and Brillouin scat- 
tering, and frequency conversion, as well as parametric 
amplification that involves trilinear boson operators |ll|: 



H — ua^a + ^ Ui ■ h\hi + g {h\h2a + h\h\i 



(6) 



The equivalence of Hamiltonians and (H) can be read- 
ily seen by applying the Schwingcr transformation 



5+ = 6l62, 5_ = hxh\, Ss = i [b\h - blb2 



(7) 



The Hamiltonians of the type of Eq. (||) are usually 
analyzed by approximating them with an exactly solv- 
able one. The solvable Hamiltonian is usually quadratic 
in the boson operators or linear in the generators of a 
classical Lie algebra. In this approach the so-called para- 
metric approximation or its variations Jl^ , [l4 | are 
often used. The basic assumption of this approximation 
is that one of the quantum modes or subsystems is pre- 
pared in a highly excited (often coherent) state, or in a 
state close to the vacuum. This approach, however, puts 
certain restrictions on the type of possible initial condi- 
tions as well as on the time span over which the quantum 
dynamics can be followed. Another type of an approach 
to the solution of Hamiltonian (||) is based on perturba- 
tion theory for nonlinear algebras p^ . This approach 
requires however the existence of a small parameter. In 
the TC model such a parameter was found in the case 
the atoms are in completely symmetric states {r ~ ^) 

The basic idea of the present work is to combine these 
two approaches through algebraic methods as applied e.g. 



|20| , pi| ) . To this end we reformulate the Hamil- 
tonian in terms of an algebra that better allows the di- 
agonalization of the Hamiltonian. This idea was already 
used by Holstein and Primakoff . They expressed the 
generators S3, S± of the su(2) algebra in terms of boson 
operators b, 5^, 



6^6, S+ 



(8) 



Here r is an index that characterizes the irrep of su(2). 
However, in [ p2| the square root in the transformation 
Eq. (||) was in the end replaced by unity, which amounts 
to applying the so-called "weak field" approximation 
^ 2r). Obviously, this approximation corresponds 
to zeroth order in the expansion of the problem with re- 
spect to parameter ^ . Transfromation (^) has also been 
applied [p3| with expansion up to second order. 

We consider here the case when the operators Aq , A± 
in Eq. (||) are generators of a polynomial deformation 
su„(2) of the Lie algebra su(2) ||2l|, |2|, |||. Numerous 
physical applications exist for polynomially deformed al- 
gebras P g l28l||, ISll M, g g M, M M. 



A particularly interesting, in view of the present prob- 
lem, application of deformed algebras was developed by 
Karassiov (see in and refs. therein). The method 
to be introduced below is an extension of Karassiov's 
method. 

We introduce here the notion of a Polynomial Alge- 
bra of Excitations (PAE). In this algebra the coefficients 
of the structure polynomials are c-numbers, rather than 
Casimir operators as is typical of polynomial deforma- 
tions. We derive an exact mapping between isomorphic 
representations of two arbitrary PAE. We also provide 
a classification of isomorphic representations of Polyno- 
mial Algebras of Excitations. In our approach classes 
of isomorphic representations are specified by the mul- 
tiplicity of the maximal and minimal eigenvalues of the 
operator Aq in the given representation. We formulate 
then an analytical approach that allows us to expand the 
Hamiltonian, when expressed in terms of PAE, as a per- 
turbation series. 

For completely symmetric states of atoms our results 
agree with those reported in Our formalism 

provides however a solution of the problem for any value 
of r, which allows us to discuss new physical effects in 
the Dicke model. 

The article is organized as follows. In the next Sec- 
tion II we discuss the irreducible representations of PAE. 
In the Section III we apply the general approach to the 
Tavis-Cummings model, and in Section IV construct the 
perturbation theory for the TC Hamiltonian and solve 
its eigenvalue problem up to third order. The generalized 
A/'-atom quantum Rabi frequency is defined for arbitrary 
quantum states of the system. In the Section V we use 
the zero order approximation for the TC Hamiltonian to 
calculate the intensity of spontaneous emission of atoms 
prepared in the state of thermal equilibrium with the res- 
onator mode. We show that the correlation of the atoms 
due to interaction with the field gives rise to the enhance- 
ment of spontaneous emission as compared to the atoms 
in the absence of resonator. In conclusion, we discuss 
possible further applications of the methods developed 
here. Technical details of the algebraic manipulations 
are given in the Appendix. 



II. IRREDUCIBLE REPRESENTATIONS OF 
THE POLYNOMIAL ALGEBRA OF 
COLLECTIVE EXCITATIONS 

The coefficients of the structure polynomial of a poly- 
nomially deformed algebra are usually expressed through 
the Casimir operators of the algebra. In this Section we 
discuss representations of a special class (PAE) of poly- 
nomially deformed algebras when the coefficients of the 
structure polynomial are c-numbers. We denote a PAE 
with a structure polynomial of order k, as il^. Formally 
Hk is an associative algebra with unity, defined by three 
generators A±,Ao. These generators satisfy two basic 
commutation relations, Eq. (|4|). As can be readily seen 



3 



from these commutation relations, [Ao,^+A_] = 0. We 
can thus assume that 

A+A^^p,{Ao)=coY[{Ao-q.,). (9) 

1=1 

Here (x) is a structure polynomial of order k, whose 
coefficients are generally complex numbers. The termi- 
nology is chosen in analogy to the structure functions 
of quantum algebras (q-deformed algebras) and the 
structure constants of the linear Lie algebras. The set 
of K real roots of the structure polynomial is denoted 
as {qi}i^i- In physical applications the operators A± of 
Eq. (H) often play the role of creation and annihilation 
operators of collective excitations. Therefore, hereafter 
the algebra Hk will be referred to as the polynomial al- 
gebra of excitations (PAE) of order k. Notice that this 
algebra is different from the algebra su„(2). All PAE's 
are completely defined by 1 c- numbers, the coefficient 
Co and the k roots {qi} of the structure polynomial. In 
the case of su„(2), however, the structure polynomial has 
some coefficients in the form of Casimir operators. 

Below we consider two elementary, but important for 
the following, examples of PAE. To begin we assume, 
without loss of generality, that cq = ±1. Indeed, in the 
case I Co I ^ 1, it is always possible to renormalize the 
generators of 11^, 

A±^\cor^A±, (10) 

such that the commutation relations Eq. (^) remain in- 
tact. Using these commutation relations it can also be 
readily seen that 

AoA± = A± {Ao ± 1) , A^A+ = (^o + 1) • (11) 

As indicated above, ^1+ and A- have the physical mean- 
ing of creation and annihilation operators of collective 
excitations (quasiparticles) , while Aq is the operator for 
the number of excitations. The most simple and impor- 
tant example of a PAE of first order, Hi, is provided by 
the well-known Heisenberg-Weil Lie algebra, viz. 

6t -> A+, 6^ h^^h^ Aq, 

00 = 1,91=0. (12) 

Here &, are the usual boson operators. For the sake of 
simplicity, in what follows we will denote the generators 
of ill by 6, 6^. The algebra ili allows us to construct the 
irrep of any other PAE of higher order k > 1 as a multiple 
tensor product of ili . 



As practical example let us describe in detail a second 
order PAE denoted in the paper by It is relevant to 
the algebra su(2) given by the commutation relations 

[S^,S±]^±S±, [5+,5_]=2^3. (13) 
It is plain that 

S+S^ = S'^-Sl + S3. (14) 

For every irrep the Casimir operator is equal to 
r [r + 1) I, where / is the identity operator. The cor- 
responding PAE of second order, §,., is constructed such 
that 

A+,S-^ A^.^Ss^ Ao; 
Co = -1, qi = -r, q2^r + 1. (15) 

Obviously, has a matrix irrep isomorphic to the irrep 
of su(2) with the same r. Notice that only in this rep- 
resentation of Sr the condition Eq.(||) is fulfilled. This 
gives us a motivation to denote the three generators of 
Sr as S^,S±. 

Below we describe a general method to construct a re- 
alization of any PAE il^ through the algebra ili . Choos- 
ing a root qj of the structure polynomial we can construct 
this realization in the form 



^0 = h]bj + g,, A+ = ^co Yl^ + q, - q,) b], 

A- = bj Co fi (b]bj + q, - q^) ■ (16) 

Here hj , 6] are the usual boson operators associated with 
the chosen root qj of the structure polynomial (x). 
The corresponding Fock vectors, i.e. the eigenvectors of 
fe^fej, are denoted as \n)j. In what follows the chosen root 
qj of the structure polynomial will be referred to as the 
pivotal root. Notice that the product in (|l^) always con- 
tains exactly k multipliers regardless of the multiplicity of 
the root. Since h^^hj has a well-defined discrete spectrum, 
the square root function is defined in the form of spec- 
tral decomposition. We specify the branch of the square 
root by choosing ^/—l = i. It is easy to see that the 
condition Eq.(||) and the commutation relations Eq.(|^) 
for generators of il^ are satisfied. 
A useful automorphism T, 

fb^ ^i-b, fb^i-b^ ^ fb^b = - {b^b + 1) , (17) 

of ill allows us to construct another realization of il^ 
through ill , 
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(-Co) Yi [1i 



q,-b]bjj, 



A_ = 



(-Co) n I'^j 



(18) 



It is worth mentioning that transformations similar to 
Eqs. (|l6|), ( [l^ ) have been introduced earher under the 
name of multiboson reahzations of Bose operators. These 
multiboson reahzations satisfy the usual boson commu- 
tation relations, [A, A'^] = 1 |||, 

Applying the realizations Eqs.(|l6|),(|l8|) in any repre- 
sentation of ill, we can construct a representation of 11^. 
In the case of realization Eq.(|l6|), an irreducible repre- 
sentation of Hk is constructed through application of the 
operator A^ to the vacuum vector |0)_,-. The finite di- 
mensional representation can be constructed in the case 
when the root Qj+i is separated from q.j by a natural 
number d. Then it can be readily seen that A%\tS)j = 0. 
Indeed, provided that A- = A^, it is not difficult to 
show that the norm of the vector j4'J_|0)j vanishes, i.e. 
j(0|AlA^[,|0)j = 0. The general construction can be 
exemplified by When r is integer or half- integer, 
d = qi — q2 = 2r + 1 is integer, and Sr has a finite di- 
mensional irrep which is isomorphic to the corresponding 
irrep of su(2). 

In the case of realization Eq.(|l8|), the corresponding 
irrep is constructed through application of the operator 
A- to the vacuum vector Using an argumentation 
similar to that given above, it can be seen that the finite 
dimensional representation can now be constructed pro- 
vided that the root <Zj-i is separated from qj by a natural 
number. 

The meaning of transformations Eqs.([l6|) andf^Sj) be- 
comes now transparent . The realization Eq.(^^ corre- 
sponds to the case when the operator A^ is a creation 
operator, while the realization Eq. ([l8|) corresponds to the 
case when A^ is an annihilation operator. 

Since the spectrum of operator b'^^bj is a set of natu- 
ral numbers and zero, the operator A^ — qj in Eq. ( p^ ) 
has a non-negative spectrum. Therefore, the argument 
of the square root function is a positive operator in the 
finite dimensional subspace, where the structure poly- 
nomial Pk;(^o) has non-negative spectrum. In this case 
the operators and A_ are Hermitian conjugated. 
In the subspace corresponding to the negative values 
of the spectrum of the structure polynomial, the argu- 
ment of the square root function in Eq.(|l^) has a neg- 
ative spectrum. The operators will be anti-conjugated, 
i.e. (A-)^ = — which is not plausible. Thus, in a 
physical problem (see (^)), we should only consider those 
irrep for which the spectrum of ^^(Ao) is nonnegative. 
The same is true for Eq. (18). 

The two relations, Eqs. (16) and (p^, between any two 
algebras il^ and ili , show that there is no principal differ- 
ence as to how exactly the meanings of the creation and 



annihilation operators of collective excitations are pre- 
scribed to the pair A±. The important point is that the 
pair exists. It is the physical problem in question which 
prescribes the meaning of operators A± and determines 
the location of the equidistant spectrum of Aq on the 
real axis. Should one be interested in the eigenvalues of 
Aq, to the right of the pivotal root qj, it is necessary to 
choose the transformation Eq.(16), while for the region 
to the left of qj , it is necessary to use the realization dlq) . 
As was explained above, we choose the region such that 
the structure polynomial is nonnegative. 

The general considerations given above can be illus- 
trated by the algebra. In the case of the realization of 

connected to the pivotal root gi = — r it is necessary 
to use Eq.(p^. This leads to the conventional Holstein- 
Primakoff representation. 



S3 = b\b, -r,S+ = V2r + 1 - b\bi bl = (^+)^(19) 
In the case of q2 = r + 1, we use the realization Eq. (08) , 



Ss^r- 462, S+ = 62 V2r -t- 1 - 6^62, 5_ = (5+)t.(20) 

The spectra of and 62^2 are limited from above by 
the value 2r, while the subspace spanned by eigenvectors 
„ (n = 0, 1, . . . 2r) of operator Aq is irreducible (cf. 
Eq. 1)). 

Obviously all the realizations of PAE constructed 
through Eqs. (p^, (18) are fully characterized by the 
dimension d of the invariant subspace and by the order 
fc_ of the left and fc+ of the right roots defining the cor- 
responding irreducible representation. Thus an irrep of 
PAE is characterized by a set of parameters d}. 
Such irrep we will denote by i?(fc_, k+,d). For instance, 
i?(l, 0, 00) means the representation of Ui, while the irrep 
of is i?(l,l,2r-h 1). 

An isomorphism between irreps of and 11'^' that 
belong to the same class i?(/c_, d), is given by 

Ao = A'o + {qj-q'y), 



A. = 



A_ = A!_ 



(21) 
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or by 



A+ = A'_ 



{q'r + - 1) - A'o: 



\ 



A. 



\ 



CO («y+9i-^o-9i) 
Co fi {^'o-1->) 



A' 



(22) 



In Eqs. (pl|),(p2|) the operator argument of the square 
root function should be taken after identical multipliers 
in the nominator and denominator are cancelled. The 
pivotal roots qj and q'^, define a vacuum vector of the 
irrep. Compare these expressions with Eqs.(p^),(p^. 

To recapitulate: in the general case two irreducible 
representations d) and R'{k'_, k'j^, d!) of il„ and 

ilJ k' are isomorphic provided that d — d' , and 



or = k'_f_^k-\- = fc' 



(23) 



For symmetric irreps, i.e. when fc+ — k-, these condi- 
tions coincide. The isomorphism is a consequence of the 
fact that, under conditions (|2^), one can choose qj{q'j) in 
such a way that the functions under the square root in 
Eqs.(^l|), (^ ) do not have zeros in the spectrum of oper- 
ator Aq, and therefore we can consider the square root as 
a single-valued analytic function. This means that trans- 
formations are invertible, analytic, and therefore define 
an isomorphism of irreps. 

The constructed transformations give us a tool to re- 
alize in terms of a simpler PAE with the same type 
of irrep. This procedure will be applied below to the 
Tavis-Cummings Hamiltonian. 



III. THE TAVIS-CUMMINGS HAMILTONIAN 
IN TERMS OF THIRD-ORDER PAE 

The interaction part of the Hamiltonian Eq. (^ can be 
expressed in terms of third-order PAE. The generators 
Mo, M± of this algebra are realized as 



(24) 



It is plain that these generators satisfy the commutation 
relations Eqs. (Q). The generators of the algebra Mq, M± 
commute with the operators 



M = a'^a + Sz+ r, S'^ = Sl + \ i^+S- H- S-S+) . (25) 

Hereafter we use the same notation M both for the 
Casimir operator and its eigenvalue, if no confusion 
arises. We show below that the eigenvalues M, r(r-f 1) of 
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FIG. 1: The structure polynomial ^3(2;) of the ] 
for r = 2. 



lM,r algebra 



the operators of Eq. (|2^) parameterize the PAE in ques- 
tion. We thus denote this PAE as M7\/^r. The structure 
polynomial of Mj\/^r can be expressed in the form 



P3 (Mo) = M+M_ 



ta(52 



02 
^3 



S3) 



Mo 



a^a (r 



- S3) ir + S3 + 1) 



M+r+2 ^ 
2 / 



The parameters of this structure polynomial are 



92 



Co = 



-1, qi 
r, qs 



M-r 



M- 



+ r + 1, 



(26) 



(27) 



and its behavior as a function of Mq is given in Fig. |^. 

We turn next to the description of finite dimensional 
irrep of M.M,r- In physical applications the parameter r 
has the meaning of collective Dicke index. This index 

runs from e(A' ) = — — to ^ with unit steps, while 
M can be any natural number including zero. Thus 53 is 
the biggest positive root of first order. If M < 2r, then 
qi > 92; if M > 2r then qi < q2] the case M = 2r 
qi — 92 corresponds to a root of second order. A typical 
plot of the structure polynomial is shown in Fig. |l| for 
these three cases. A finite dimensional representation of 
Mjv/.r (where (M_|_)'l' = M_) corresponds to the positive 
spectrum of (Mq). The spectrum is limited from the 
right by ^3 and from the left by qi or 52- Notice that 
the number M is conserved. Therefore this number is 
determined by the initial state. The different values of 
M and r define different algebras Mjvf.r, whose single 
physical finite dimensional representation we will call a 
zone. The case M < 2r corresponds to nearby zones. 
The two largest roots are qi and 53 and the irrep has 
the type i?(l,l,M-|- 1). Consequently, the well-known 
weak- field limit corresponds to nearby zones. 

The case M > 2r corresponds to remote zones. The 
two largest roots are q2 and (73, and the corresponding 
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we obtain 




FIG. 2: Approximation of the structure polynomial P3{x) 
by tlie parabolic P2{x) for r = 2. a) M = 1; b) = 4; c) 
M = 12. 



irrep is of the type i?(l, 1, 2r + 1). Notice that the region 
M >> 2r is usually called the strong-field limit. 

In the special case 2r — M, referred to intermediate 
zone, the algebra M.M,r possesses an irrep of the type 
R{2, l,2r + 1). It is the only irreducible representation 
that principally differs from all the others. 

As indicated above, the simplest PAE with irrep of the 
type 1, d) is Sf (we use here f to distinguish it from 
the (physical) collective index r) . It would be convenient 
to solve the eigenvalue problem in terms of the simplest 
algebra Sf. Notice that in any finite dimensional irrep 
of PAE characterized by R{1, 1, d), the structure polyno- 
mial of the algebra can be approximated by a parabolic 
curve. This is shown in Fig. ^ (for M 7^ 2r) for the 
structure polynomial of MM,r- The larger (smaller) M is 
in comparison with 2r, the better is the approximation. 
However, for AI w 2r the approximation is not satisfac- 
tory (see Fig. ^). In the regions where the approxima- 
tion is adequate, it is then not difficult to diagonalize 
the operator V = S+ + defined in terms of genera- 
tors of the conventional su(2) algebra. The latter has a 
parabolic structure polynomial. In Fig. |2| the two roots 
of p2 {x) are chosen to be equal to the corresponding two 
roots of P3{x). The choice of cq in P2{x) will be explained 
below. Thus in the cases of nearby or remote zones, it is 
convenient to study the problem in terms of algebra Sf. 
The approximation illustrated in Fig. |^ indicates that the 
TC problem can be solved via an appropriate perturba- 
tion theory. 

To begin with we consider the transformation of Ma/^^ 
to Sf for the case of remote zones. The dimension of a 
remote zone is 2r + l, and the algebra §,~ should be char- 
acterized by f = r. The finite dimensional irrep of §,~ 
is isomorphic to the corresponding irreducible represen- 
tation of the atomic subsystem. For the pivotal root qj 
we choose the largest root that bounds the irrep of MM,r 
from the right (the in ( p7^)) , while as the root q'^, we 
take the qi of §r from Eq. (Ea). Applying the mapping 



Mn 



M-r 



M - r + 1 - 53 , 



M_ = ,/(Af-r + l-53 5+. 



(28) 



The spectrum {m} of the operator 53 belongs to the re- 
gion —r < rh < r, conse quen tly the argument of the 
square root function in Eq.(pq) is positive in the remote 
zones (M — r > r). The relations Eqs.(p^) express the 
generators of algebra M.M,r as analytic function of the 
generators of the algebra, they thus allow us to ap- 
proximate the more complex algebra MM.r of third order 
by a simpler algebra of second order. 

Before we begin to study this approximation we con- 
struct a realization of in terms of boson and spin vari- 
ables. iFrom Eqs.(|2|),(|5l), and <^ it follows that in 
remote zones 



Si — S3, S+ — 



1 



a ■ 5+ , S- = ( S- 



t 



(29) 



Notice that the subspaces that correspond to remote 
zones do not contain the vacuum state of the field. It 
is also worth mentioning that the matrix representation 
=a is 6n,n+i in any remote zone. This 



of operator , 

operator has been considered before as a earlier in terms 
of phase operator pj| . 

We turn now to the nearby zones M < 2r. For this 
region the mapping of the algebra MM,r to the algebra 
Sf is realized through procedure similar to that described 
above for remote zones. Notice that the dimension of 
nearby zones is d = qs — qi = M + 1, and therefore 
f = ^ (cf. Fig. ^). For the nearby zones there is no 
simple correspondence to a physical representation like 
that of Eq. ( p9|) . The particular choice of the pivotal root 
q'j, is unimportant and we use the same choice as before. 



Applying Eq.(|22D we obtain 



S3, 



S^ 



4r-M 



1-S. 



M. 



(30) 



Since all the eigenvalues of the operator ^3 belong to the 
interval — f to f, the argument of the square root function 
docs not have zero eigenvalues in the nearby zones. The 
realization of through spin and boson variables is then 
given by 



M 
~2 



S. 



x/r + 1 - 53 
1 



S+ a, 



S- 



+ 1 - ^3 



(31) 



Notice that the nearby zones do not contain the eigen- 
vector |r, r) of S3. 
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To clarify the structure of intermediate zone we choose 
as the pivotal root as the third root of ^3(0;) of Eq.(27), 
and apply the transformation {UM. We thus obtain 



Mo = 93 - bibs -1 = M±L- bibs, 



M. 



M + 1- bibs I ( 2r + 1 - blk 



M +1 



bibs) (2r+ 1-463) bl (32) 



Here 63, 63 are generators of the algebra ili 



In the region M > 2r the multiplier (2r + 1 — 6363 ) m 
Eq.(^2[) vanishes first. This corresponds to remote zones. 
In the region M < 2r the multiplier that vanishes first 
is (M + 1 — 6363). In the intermediate zone M = 2r we 
have 

M+ 63 (2r + 1 - 6;^63), = (2r + 1 - bibs) bl (33) 

In the intermediate zone we thus obtain a special real- 
ization of the TC Hamiltonian, 



H{M=2r}^i2r+^)ibs+bl)- 

blbs{bs + bl) + ibs + bl)blbs 



(34) 



The domain of the quantum space for this Hamiltonian 
is specified by the condition ns < 2r. 



IV. DIAGONALIZATION OF THE 
TAVIS-CUMMINGS HAMILTONIAN 

Let us introduce an AA-atom generalization for arbi- 
trary values r and M of the well-known quantum Rabi 
frequency such that 



2jM-r + i M>2r 



(35) 



4r-M+l 



, M <2r 



For r = =^ our definition agrees with that used in|l^, 

|r^ . Introducing a small parameter a = (^$1^) ^ we 
can rewrite the realizations of AI± in nearby and remote 
zones (see Eqs.(|2^),(^o|)) in the form 



M4 




1 - a ( 53 - ^ ) , M- = M\ 



(36) 



The diagonalization problem for the operator Eq. (|^) can 
now be solved in each zone by means of perturbation 
theory with respect to the small parameter a. One can 
show that the eigenvalues of the argument of the square 
root function in Eq. ( ^ ) are less then unity. Hence we can 
expand the square with respect to a, and find thereby for 
the interaction part of the Hamiltonian 



J 



V 



(5+ + 5_)-f ((^3-i)5+ + ^_ (^3-i)) - J,(f)^((^3-i)'5+ + 5_ (53 -i 



E 

n— 3 



(2n-3)!! j aiSs-i) 



S+ + S. E 

n— 3 



(2«-3)!! / ajSs-^) 



(37) 



In the interaction representation the Hamiltonian coin- 
cides with V, we only need to diagonalize the latter. Up 
to third order in a we find that 



V 



(3) 



where the V^""^ are terms of n-th order in a, and are given 
in the Appendix |A| In the Appendix we also show that 
unitary transformations Uk, k = 0,1,2,3, which bring 
the interaction operator into diagonal form: 



v = u-v-u-^^nR-Ss<i + ^ 



5Sl 



3f {f 



1) + 1 



U = UsU2UiU„. 



(38) 



The spectrum of the operator V as given by Eq.(|3q) 
agrees with the results of |l^, ^ for the symmetric states 
of the atoms. 



We compared the third order solution Eq.(|38|) with 
the exact numerical diagonalization of V and found that 
the result ( |38| ) is very accurate, especially for increasing 
values of \M — 2r\. The results of this comparison are 
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FIG. 3: Deviation of the eigenvalues of V from their numer- FIG. 5: 
ical values in zero (solid line) and second (dashed line) order y for r 
in a, for r = 6. 



The dependence on M of the maximal eigenvalue of 
= 3, as calculated numerically and from Eq. 



o Exact numerical solution 

X Algebraic approach 

• Tavis-Cummings's solution 

2500-1 



100- ••• 



26 30 36 j 40 46 




FIG. 4: Energy levels of V in the ascending order for r = 25, 
as calculated numerically (circles), from Eq. (38|) (crosses), 
and as is given in M (squares), a) M = 50; b) M = 2525. 



Finally we consider the A/'-atom dressed states. In 
other words we introduce a representation in which the 
zero order Hamiltonian is diagonal. This representation 
is given by transformation Uq. For J\f — I all the higher 
order terms in Eq.(|37|) (higher than zero order) vanish, 
while the eigenstates of the zero order Hamiltonian coin- 
cide with the dressed states Q of the Jaynes-Cummings 
model. We can thus call the eigenstates of the zero or- 
der Hamiltonian the dressed states of the A/'-atom model. 
Notice also that in the remote zones the Rabi frequency 
Eq.(^5|) does not depend on the Dicke index r. If we con- 
sider only the zero order terms of _ff , it is convenient to 
combine all the remote zones that have the same Rabi 
frequency into a remote superzone, whose dimension is 
2^. Introducing an operator for the total number of 
quanta in the atom-field system (cf. the definition of M 
in Eq.(|2F 



t o 

K = aTa + S3 + —, 



(39) 



shown in Fig. |^. 

In their original article [|| Tavis and Cummings also 
found an approximative analytical expression for the 
spectrum of the interaction operator V. We compare 
the results of Ref. with our analytical and numerical 
solutions in Fig. ^. It is evident that the Tavis and Cum- 
mings solution is only accurate for very large and very 
small values of index j. 

The Fig. I compares the energies calculated numeri- 
cally and in accordance with the analytical solution Eq. 
(^). In the intermediate region of M the curves for 
nearby and remote zones overlap and coincide thus pro- 
viding still satisfactory correspondence to the exact so- 
lution. However evidently the expansion for the remote 
zone breaks down in the nearby zone and vice versa. This 
means that the classification of zones introduced in this 
paper is indeed adequate. 



we can define the remote superzone as follows. A re- 
mote superzone contains all eigenvectors of K that have 
the same eigenvalue, provided it is larger than Af. In 
Fig. ^ we show the definition of the zones in the Tavis- 
Cummings model for A/" = 4. 

It can be readily seen that, in the Hilbert subspace cor- 
responding to the remote superzone, the unitary transfor- 
mation Uq can be factorized into a product of J\f single- 
particle transformations that are the dressing operators 
for the Jaynes-Cummings problem (cf. Eq.(p9[)), 



Uo 



J|exp 



1 



VaJo+T 



1 7 t 



1 



(40) 



^From the discussion following Eq.(p9|), it follows that 
in each remote superzone the field operators 



Y^oto+l 
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r-2 
6(2) = 1 
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0(1) = 2 0(0) =1 



the system. 

In this section we study a contribution of cooperative 
effects into the rate of spontaneous emission generated by 
the atom-field system. The atom-field system is assumed 
to be prepared in the state of thermal equilibrium. This 
state is described by the canonical Gibbs ensemble with 
the thermostat temperature T. This means that under 
the " system" we imply JV atoms strongly coupled to the 
resonator mode. Under the "thermostat" we imply the 
surrounding environment, for instance, cavity walls taken 
at the temperature T. In the state of thermal equilibrium 
the exact atom-field density matrix should be defined 
using the Hamiltonian Eq.(|l|). However, for the exact 
density matrix the analytical analysis of the transition 
probability, if possible at all, would be highly technically 
involved. 

We demonstrate that nontrivial physical results for the 
intensity of spontaneous emission of two-level atoms 
placed inside the cavity can be already obtained for 
the zero order approximation of the exact Hamiltonian 
Eg. (pT\) . The thermal state is given by 



Pth = cxp 



Ho 



S+ + S- 



kT 



(41) 



FIG. 6: Zones structure for A/" = 4. iS'^-* is the spectrum of 
the zero order Hamiltonian. K is an eigenvalue of the corre- 
sponding operator (^), and G(r) is a number of equivalent 
irreps of su(2) defined below. 



and 



commute. This means that the system 



behaves almost like a semiclassical one. 



V. ENHANCEMENT OF SPONTANEOUS 
EMISSION IN THE RESONATOR DUE TO 
COLLECTIVE EFFECTS 

In the previous section we developed an algebraic ap- 
proach to the Tavis-Cummings model. We introduced 
the operators S± describing collective excitations in the 
atom-field system. In terms of these operators we con- 
structed a perturbation series for the Tavis-Cummings 
Hamiltonian Eq.(|^). The derived perturbation series 
gives us a tool to distinguish and classify cooperative 
(multiparticle) effects of different orders that are involved 
in calculations of different physical observables charac- 
terizing the atom-field system. In the preceding sec- 
tion we constructed a Hamiltonian Eq.(p7|) correspond- 
ing to the zero order approximation for the Hamiltonian 
Eq.(|l]). This simplified operator depends on multiparti- 
cle Rabi frequency Eq. ( p5| ) , which depends on the number 
of atoms in the cavity. Therefore the simplified Hamilto- 
nian Eq. (|37]) allows to account for cooperative effects in 



Here Z is a normalization factor. 

We show here that the intensity of spontaneous emis- 
sion of the system comprised of Af atoms strongly coupled 
to the resonator mode and prepared in the thermal state, 
can be greatly enhanced at a certain temperature. This 
amplification results from high correlations in the atomic 
subsystem. Similar effect exists for Dicke's superradiant 
state iQ. This state is prepared by a short laser pulse. 
Therefore the setting of the Dicke's theory is quite dif- 
ferent from our considerations of the stationary state of 
thermal equilibrium. 

When calculating the rate of spontaneous emission (or 
the intensity proportional to this quantity) we merely 
follow the ideas of Dicke's paper (see e.g. Ql). According 
to this theory, the rate of spontaneous emission in the 
system is proportional to the average of the square of the 
atomic dipolc, viz. 



/ = /o • {S+S-) =h-Tr {pth ■ S+S-} . (42) 

It is convenient to calculate the average in the dressed 
states basis, where 



1 ~ ~_i 
Pth —UopthUg = exp 



Hp + ^rSs 
kT 



(43) 



Taking into account Eq.(|T^ along with the fact that, for 
an arbitrary zone. 



53 = (f-r) + 53, 



(44) 
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it can be shown that in the dressed states basis 
f/o(5+5-)C/o-i=f(2r-|f + i)- 
(r + 1 - r}{S+S-) -1{SI- 2Si + Si) . (45) 

FinaUy, the intensity of spontaneous emission is given by 



AT 

oo — 



M=Or=e 



e ET [f (2r - |r + i) + \m^] . (46) 



Here e = ^—^4^, and G(r) = is the 

number of equivalent representations with the same r. 

Let us consider the intensity per atom, i.e. Ii = I /J\f. 
This intensity consists of two terms, i.e. the first is given 
by a single-particle contribution Isingie and the second 
one proportional to the two-particle correlation function 
Cor, /i = Isingie + h{J^ ~ l)Cor. They are found to be 

Isingie = hjj{Y.^+^-) = {\ + jj^^'i)) ' 

i 

It is plain that in the absence of the cavity, the correlation 
function vanishes and the only contribution to I\ is given 
by the first term. 



(48) 



Isingie = /c/ = /o I 1 + 6 



The contribution of any remote superzone can easily be 
found due to the factorization property of the operator 
Eg. (pO[) . If we denote the trace over states that belong 
to the same remote superzone by subscript if, we obtain 



(2cosh(il-))^-^(sinh(il-))^ (49) 

In Fig. 1^ we compare the spontaneous emission Ii in the 
presence of a cavity against the intensity Id of M atoms 
in the absence of cavity. 

Notice that if the number of atoms is big enough, the 
intensity of radiation exhibits a high maximum. In a 
cavity at low T, the cluster of M two-level atoms emits 
much more intensively than it does in the free space. It 
should be possible to drive the system to thermal equilib- 
rium at the temperature where the spontaneous emission 
exhibits maximum. The marked amplification of sponta- 
neous emission should be observed in cavity experiments. 

Concluding this section, we recapitulate our main re- 
sults. We consider spontaneous emission of the system 
comprised of N two-level atoms strongly coupled to the 
cavity mode and prepared in the state of thermal equi- 
librium. In the absence of the cavity the atoms in the 
thermal equilibrium would be uncorrelated. In this case 



2.0 



1.5 



1.0 



0.5 



0.0 



(1) 
(2) 
(3) 
(4) 



10 20 30 40 50 60 70 

kT 



FIG. 7: The intensity of spontaneous emission per atom 
(in units of Jo) versus cavity temperature. The curve (1) 
is the classical result given by Eq.(^8[). The curves (2,3,4) 
correspond to M = IG, M = 50, M = 100, respectively, and 
^ = 10. 



the spontaneous emission would be described by the con- 
ventional formula Eq.(48). For high-Q resonators the 
strong coupling to the resonator mode should necessarily 
be taken into account. We demonstrate that in this case 
the intensity of spontaneous emission can be be greatly 
enhanced. This phenomenon can be explained by addi- 
tional correlation between atoms established by the cav- 
ity mode. To analyze the effect analytically we have re- 
placed the exact Tavis-Cummings Hamiltonian Eq.(|l|) by 
its zero order approximation derived in the previous sec- 
tions. This allowed us to represent the intensity of spon- 
taneous emission in the simple analytical form Eq. (^6[) . 
It is appropriate to emphasize once again that the zero 
order approximation of the Hamiltonian contains strong 
coupHng and, thus, describes cooperative effects in the 
atomic subsystem. This is the consequence of the fact 
that the operators S± describe collective excitations in 
the atom-field system. 



VI. CONCLUSION 

In this work we solved the Tavis-Cummings problem 
by applying the technique of polynomially deformed alge- 
bras. We constructed the transformations that map one 
polynomial algebra of operators onto another. This al- 
lowed us to reformulate the problem in terms of a simpler 
algebra of second order, §,., and develop a specific pertur- 
bation theory. Our results have a significant advantage 
over the so-called linearization approximation, i.e. the 
case when the Hamiltonian is linearized in terms of the al- 
gebra ill ■ In this latter approach a structure polynomial 
of a higher order is approximated by a polynomial of first 
order. This method allows only to calculate the lowest 
or highest eigenvalues, and the corresponding eigenstates. 
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It does not take into account the finite dimensionality of 
the representation. The parabohc approximation devel- 
oped in this work provides in this respect a significant 
advantage because it allows to construct a finite dimen- 
sional representation for the problem. We were able to 
find analytical expressions for all the eigenvalues of the 
Hamiltonian up to third order in the small parameter 
a. For the nearby zones we showed explicitly how the 
collective quantum Rabi frequency depends on Dicke in- 
dex r. Since this index characterizes the symmetry of 
atomic states, the result has significant physical implica- 
tions. The dependence on atomic symmetry is revealed 
already in zeroth order in the perturbation expansion. 
Employing our methods we found an interesting new ef- 
fect, amplification of spontaneous emission of thermal A/"- 
atom states due to collective effects. We expect that this 
phenomenon can be observed in cavity experiments. It 
is worth mentioning that the applicability of the method 
developed in this paper can be extended to many other 
problems, including Bose-Einstein condensation, multi- 
photon interactions in the micromaser, and multimode 
interaction of an electromagnetic field with matter. Our 
methods allow in particular to address the problem of 
collective and dressed states in the mentioned physical 
systems. 



V 



(1) 



V 



(2) 



where 



and 



~ 5+ + S- ~ 
^3 7: ^3 



V(3) = - 



1 fa\^ 

2 V 27 



B = S^Sx + SxSg 



S+ + S- 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



1. The zero order transformation Uo 
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As known from the theory of su(2) algebra, the oper- 
ator Sx can be diagonalized by the transformation 



Uq = exp 



S+-S- 



= exp 



Employing this transformation we obtain 



(A8) 



APPENDIX A: SIMILARITY 
TRANSFORMATIONS 

We look for similarity transformations Uk that diago- 
nalize the Hamiltonian in different orders of a. 



V^n^ = 2^3 ( 1 



1 /a\2 



2 V4 



^(2) 



^(3) 



SxSsSx 



1 /'a\3 r~ 



(A9) 



(Al) 



where fc = 0, 1, 2, 3. In Eq. ( Al)) only the terms of order 
n in the small parameter a are present. For fc = 0, the 

(n) 

term V^i should be replaced by the corresponding term 



in Eq. 



Regrouping the terms we obtain 



-(-)' 

2 V4/ 



(A2) 



2. The first order transformation Ui 

It can readily be seen that the transformation 

Ui = exp [ax ■ Dl] , Dl = -i \S3Sy + SySsl (AlO) 

diagonalizes the operators in the first order. In the diag- 
onalization one needs the commutators 
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S3,D1 =B, 



S3,D1 



,D1 



45*3 - 2Sl 



,D1 



,D1 



1 

3-4 
S3B + BS3 - S3 



S3Sj3,Dl\ = 2Sl - 2sf) + (^3^.)' + (5.53)' • 



(All) 



Then up to third order in a, 



-2ax\l 



1 /a\2 



1 /a 

2 V4 



4 (aa;)^ 5*3 ( 5*^ - 2~Sl - ^ 



(ax) 



3 



{ (452 - 1) B - 8 j} (A12) 



and 



iJ.v'i^ (C>i)" = f (l + Hf)')•^- 
(52-25|-i)+a3(|)^ 

{ (452 - 1) B - 8 j} , 



a^a; . 253 ( S"" - 2S? 



where J = S3B + BS3 - 838x83. 



(A13) 



^(2) 
Vo 



a' (f )' • {5.^3 (28^ - 481 - + (2^2 - 45| - S^) ^s^,} 



In third order the operator Vq remains unchanged after and find then that 
the transformation, i.e. V^^^ — Vq^\ In order to calcu- 
late Vq up to third order, we take into account that 



_lil+(fQ_^l / ...2X 



vr = vr - 4 



52 - 25? 



Vf ^ = V^^^ -Hj)'B + Ulf{ (4^^ - 1) • B - 8J} + 
(f )' • {^.^3 (252 - 45| - (5,)') + (^252 - 45| - (5,)') 838,^ 



(AM) 



(A15) 



3. The second and the third orders transformations 



To find the second order transformation we rewrite ^ ~ ~ 2 (2^) "^^"^^ '^3-^x + 83 (^8^ — 8^ — 1 



(2) 

Vi in a symmetrized form. 



(A16) 
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where Lx 
then given by 



si+st 



The diagonaUzing transformation is 



U2 = exp <; i- 



1 fa 



(A17) 



si -si 



where Ly = 
we obtain (see ( A13 )) 

a.vS«>(S.)-' = vr + i(f) 



. Keeping the terms up to third order 



(.A18) 



The transformation (A17) does not change the expres- 



/ 2) ('3') 

sions, given above iorV^ ' and V^/ ' , and we find that 



(3) 

are no diagonal terms in , which would contribute to 
the spectrum of the Hamiltonian, we do not give here the 
fairly complicated form of operator O. The final diagonal 
form for the interaction V is thus given by 



V = nS3<l 



bSi ~ 3r {f + 1) + 1 }. (A20) 



To recapitulate, we introduced four transformations 
Uk, fc = 0, 1, 2, 3, which successively diagonalize the in- 
teraction operator in the Tavis-Cummings Hamiltonian 
up to third order with respect to the small parameter 



with f2fl the generalized Rabi frequency 



of Eq. (|35|) , such that 



(A19) 



(3) 

Diagonalization of V2 can be performed in a similar 



way with an operator U3 = exp 



■O 



As there 



V = UVU^ 



U = U3U2U1U0. 



(A21) 
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